Theory of two-dimensional macroscopic quantum tunneling in YBa 2 Cu307_«5 
Josephson junctions coupled to an LC circuit 



Shiro Kawabata, 1,2 ' 3 Thilo Bauch, 2 and Takeo Kato 4 
1 Nanotechnology Research Institute (NRI), National Institute of Advanced 
Industrial Science and Technology (AIST), Tsukuba, Ibaraki, 305-8568, Japan 
2 Department of Microelectronics and Nanoscience (MC2), 
Chalmers University of Technology, S-41296 Goteborg, Sweden 
3 CREST, Japan Science and Technology Corporation (JST), Kawaguchi, Saitama 332-0012, Japan 
4 The Institute for Solid State Physics (ISSP), University of Tokyo, Kashiwa, Chiba, 277-8581, Japan 

(Dated: November 4, 2009) 

We investigate classical thermal activation (TA) and macroscopic quantum tunneling (MQT) 
for a YBa2Cu307_<s (YBCO) Josephson junction coupled to an LC circuit theoretically. Due to the 
coupling between the junction and the LC circuit, the macroscopic phase dynamics can be described 
as the escape process of a fictitious particle with an anisotropic mass moving in a two-dimensional 
potential. We analytically calculate the escape rate including both the TA and MQT regime by 
taking into account the peculiar dynamical nature of the system. In addtion to large suppression 
of the MQT rate at zero temperature, we study details of the temperature dependece of the escape 
rate across a crossover region. These results are in an excellent agreement with recent experimental 
data for the MQT and TA rate in a YBCO biepitaxial Josephson junction. Therefore the coupling 
to the LC circuit is essential in understanding the macroscopic quantum dynamics and the qubit 
operation based on the YBCO biepitaxial Josephson junctions. 

PACS numbers: 74.50.+r, 03.67.Lx, 74.72.-h, 85.25.Cp 



I. INTRODUCTION 

Macroscopic quantum tunneling (MQT) has become a 
focus of interest in physics and chemistry because it can 
provide a signature of quantum behavior in a macroscopic 
systemic Among several works on MQT, Josephson 
junctions have been intensively studiedi^il Heretofore ex- 
perimental tests of MQT were focused on low-T c super- 
conductor Josephson junctions. 

Renewed interest in MQT occurred after the recent 
experimental observations of MQT^M i 10 ' 11 and en- 
ergy level quantization (ELQ)£ i 12 i 13 i 14 in high-T c su- 
perconductor Josephson junctions, e.g., YBa2Cu307_a 
(YBCO) grain-boundary biepitaxial junctions [see Fig. 1 
(a)p£ and Bi 2 Sr 2 CaCu208+5 intrinsic junctions^ High- 
Tc Josephson junctions, characterized by a high Joseph- 
son plasma frequency uj p (up to several THz), exhibit a 
crossover from the thermal activation (TA) to MQT at 
relatively high temperatures in comparison with low-T c 
systems. 

These intriguing findings have definitively opened up 
the way to quantum systems based on the <i-wave sym- 
metry of high-T c superconductors. The d-wave or- 
der parameter can be used to create naturally degen- 
erate two-level systems which can offer significant ad- 
vantages for quantum computation i 17 i 18 i 19 i 20 i 21 i 22 How- 
ever, one of the main arguments against high-T c qubits 
based on systems with d-wave order parameter symme- 
try was the presence of low energy excitations, i.e., nodal 
quasiparticles and zero-energy Andreev bound states 
(ZESs), destroying quantum coherence. Recent theoret- 
ical works suggest that nodal quasiparticles and ZESs, 
respectively, give super-Ohmic and Ohmic dissipation on 



MQT , 23 - 24 - 25 - 26 - 27 - 28 - 29 

The YBCO biepitaxial grain-boundary junctions 
which were used in MQT and ELQ experiments^^ have 
quite novel structure and several fascinating advantages 
for the coherent property over other types of high-T c 
junctions. First, the relative orientation of the d-wave 
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FIG. 1: (Color online) (a) Schematic of a biepitaxial YBCO 
junction and (b) the extended circuit model for (a), including 
the bias current 7 cxt , the stray capacitance Cs and kinetic 
inductance Ls- The Josephson junction with the critical cur- 
rent Ic and the capacitance Cj is formed at the boundary 
between (001) and (103) YBCO electrodes. <f> and (j>s are the 
phase difference across the Josephson junction and the stray 
capacitance Cs, respectively. 
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order parameter between two YBCO electrodes can be 
designed artificially ! 30 ' 31 Therefore it is possible to make 
a junction in which a lobe of the d-wave order parameter 
of one electrodes is facing a node in the other electrodes 
as shown in Fig. 1(a). In such configuration, super- 
Ohmic dissipation resulting from the nodal quasiparticles 
can be suppressed] 23 ' 26 ! 28 Moreover the junction has an 
additional tilting of one electrode [(103) YBCO] with re- 
spect to (001) YBCO electrodes [sec Fig. 1(a)]. In such 
tilt-junctions, the formation of ZESs causing Ohmic dis- 
sipation can be inhibited as was predicted by Golubov 
and Tafuri^ 2 . 

Recently, it was found from ELQ experiment that 
the measured bias-current I cxt dependence of the reso- 
nant frequency deviates significantly from what one ex- 
pects from the Josephson plasma frequency ui p in a con- 
ventional single Josephson junction.— The resonant fre- 
quency is a factor of five less than what one would expect 
from the estimated values of the Josephson critical cur- 
rent Ic and the junction capacitance Cj. This fact could 
be explained by taking into account the existence of a 
large kinetic inductance and stray capacitance coupled 
to the junction^ 

In the biepitaxial junctions, the film of the (103) 
YBCO electrode is oriented in such a way that the elec- 
tric transport has a large component in the c-axis di- 
rection [see Fig. 1(a)]. So the kinetic inductance Lg of 
the (103) electrode can be significant and is much larger 
than the Josephson inductance of the junction Lj*£ An- 
other important effect which was observed in ELQ ex- 
periments is the influence of the stray capacitance Cs of 
the SrTiOa (STO) substrate with a huge dielectric con- 
stant [see Fig. 1(a)]. Hence a large part of the measured 
capacitance is not due to the junction interface but to a 
distributed stray capacitance Cs which in effect shunts 
the junction^ 

The influence of and C's on the macroscopic dy- 
namics is inevitable in the biepitaxial junctions and can 
be taken into account by an extended circuit model [see 
Fig.l(b)] In Fig. 1 (b), 4> is the phase difference 

across the Josephson junction and (f>s = (2ir /<& )IsLs+(f> 
is the phase difference across the capacitor Cs, where 
Is is the current through the inductor Ls- As will be 
mentioned later, the addition of the LC circuit results 
in a two-dimensional (2D) potential U((f>, 4>s) an d an 
anisotropic mass which make the dynamics much more 
complex than for an ordinary single junction. In Rcfs. 
fl^l33ll3"4j we have theoretically investigated the classical 
resonant activation based on the extended circuit model, 
and found that the bias-current J ex t dependence of the 
resonant frequencies can be quantitatively explained by 
the normal modes (upper and lower resonant mode) of 
the system. Those normal modes, under certain condi- 
tions, can be completely different from the bare Joseph- 
son plasma frequency tu p . Therefore the extended circuit 
model well describes the ELQ experiment in the YBCO 
biepitaxial junctions. 

In contrast to ELQ experiment, it was experimen- 



tally confirmed that the thermal escape process above 
the crossover temperature T co can be quantitatively ex- 
plained by the single junction model without the LC 
circuit, i.e., one-dimensional (ID) models Therefore a 
natural question to ask is can the extended circuit model 
explain both the TA and MQT escape dynamics in the 
YBCO biepitaxial junction consistently? In this paper, in 
order to answer the important open question, we explore 
the validity of the extended circuit model for both the 
TA and MQT escape processes by extending the previous 
classical theory! 33 ' 34 Then we will show that the presence 
of the LC circuit has negligible influence on the TA escape 
rate, but the MQT escape rate is suppressed considerably 
due to the coupling to the LC circuit. This behavior 
is nicely consistent with a recent experimental result of 
the temperature dependence of the TA and MQT escape 
rate.— Therefore, the presence of the LC circuit is funda- 
mental and essential in understanding the macroscopic 
phase dynamics in the YBCO biepitaxial junctions. 

This paper is organized as follows. In Sec. II, we show 
the Lagrangian describing the YBCO Josephson junction 
coupled to the LC circuit. We also discuss the anisotropy 
of the mass and the two-dimensional potential profile for 
this model. In Sec. Ill, we derive an effective action 
from the Lagrangian and show that the system can be 
mapped into a one-dimensional model. The crossover 
temperature T co , the TA, and the MQT escape rate are 
calculated in Sees. IV, V, and VI, respectively. We try 
to compare theoretical results with recent experimental 
data of the YBCO biepitaxial Josephson junction in Sec. 
VII. Finally, we summarize our results and draw future 
directions in Sec. VIII. 



II. MODEL 

In this section we derive the Lagrangian for the Joseph- 
son junction coupled to the LC circuit and discuss the 
anisotropy of the mass and two-dimensional potential 
structure of this model. The Hamiltonian of the circuit 
[Fig. 1(b)] can be written as 
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where Qj = Cj(^ /2n)(d^/dt) and Q s = C s ($ q /2-k) 
(d<f>s/dt) are the charge on the junction and the stray 
capacitor Cs, respectively, Ej = (h/2e)Ic is the Joseph- 
son coupling energy, and $0 = h/2e is the flux quantum. 
The Lagrangian is then given by 
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where 7 = I cyLt /Ic and 77 = 2TrI c L s /§o = L S /L J0 with 
A/o = $o/2nIc being the zero-bias Josephson induc- 
tance. The Lagrangian C describes the quantum dy- 
namics of a fictitious particle with an anisotropic mass 
moving in a two-dimensional tilted washboard potential 
U{<t>i4>s)- Therefore the escape paths traverse a two- 
dimensional landscape, so the macroscopic dynamics in 
this model becomes more complicated than that in the 
simple one-dimensional model. 

The two-dimensional potential profile U((f>,4>s) is 
shown in Fig. 2. The mean slope along the (f>s di- 
rection is proportional to the bias current 7 cx t, and the 
mean curvature perpendicular to the diagonal direction 
(4> = (f>s) is due to the inductive coupling between the 
Josephson junction and capacitance Cs characterized by 
77. The local minimum point (</> m ,0™) and the saddle 
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FIG. 2: (Color online) The two-dimensional potential profile 
U((f>, <f>s) of a YBCO Josephson junction coupled to a LC cir- 
cuit for 7 = Icxt/Ic = 0.92 and 77 = Ls/Ljq — 7. Initially 
the junction oscillate around the local minima at (<f> m ,({fg). 
If 7 is increased, the junction eventually switches to the fi- 
nite voltage state by escaping either thermally or via MQF 
through the saddle point at (</>*, </>g). 



5^*, 0$) is given from Eq. ([3]) as 

b m ,<t%) = (sin" 1 7,777 + sin" 1 7), 
0s) = (7r-sin -1 7,7r-sin -1 7- 



777) 



(4) 
(5) 



Then the potential barrier height V^ 2D is given by 



v 2D = u. 



= Ej (27 sin -1 7 - 7T7 + 2^/1 - 7 2 ) , (6) 

and is a decreasing function of 7. Importantly, the poten- 
tial barrier height Vq 2D does not depend on the LC circuit 
parameters Ls and Cs and the expression of V^ 20 is iden- 
tical with the barrier height for usual single Josephson 
junctions, i.e., one-dimensional models 

On the other hand, the overall potential shape can be 
controlled drastically by changing L5. In the case of 
small Ls or equivalently large coupling strength 77" , the 
side well becomes shallower and the potential is almost 
confined to the diagonal direction (0 = <fis)- On the other 
hand, when the coupling strength 77 _1 is decreased, corre- 
sponding to a large Ls, the confinement to the diagonal 
direction tends to be weak. This makes the escape from 
a metastable well much easier. In Sees. V and VI, we 
quantitatively investigate the Cs and Ls dependence of 
the thermal and quantum escape rate. 

Next we derive an approximate expression of U(<f>, 4>s) 
f° r 7^1- By introducing the new coordinate (x,y) = 
(4> — (j) m , <ps — 0™) and assuming 7 < 1 , we can rewrite 
the Lagrangian as 
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In 



this equation, Um(x) = Ej — x 3 /6 + y{T— j)/2x 
(l/2)Mujp(x 2 — x 3 /xi) is the potential of the Josephson 
junction without the LC circuit, where x\ — 3 \J 1 — 7 2 is 
the second zero of U\v>(x) and uj p = uj p o{1 — 7 2 ) 1 / 4 is the 
Josephson plasma frequency with oj p q = ^J2ttIc/WqCj = 
\j \JLjqCj being the zero-bias plasma frequency. 



III. DERIVATION OF EFFECTIVE ACTION: 
MAPPING TO A ONE-DIMENSIONAL MODEL 

Here we derive an effective action from the Lagrangian 
Eq. (J7J) of the extended circuit model. Note that, in the 
ELQ experiment, the quality factor Q of the junction 
was found to be sufficiently large (Q m 40)^2 so we can 
safely neglect the damping effect resulting from external 
environments other than the LC circuit. 
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By using the functional integral method^^ the parti- 
tion function Z of the system can be written as 



Z = 



J Vx(t) J Py(r)exp < 



dT£[x,y] , (9) 



where the functional integral is performed over all the 
periodic paths with the period h(3. In this equation, 
C[x,y] = (M/2)x 2 + (m/2)y 2 + U(x,y) is the Euclidean 
Lagrangian and [3 = 1/ksT. The Lagrangian is a 
quadratic function of y and the coupling term between 
x and y is linear, so the functional integral over variable 
y can be performed explicitly by use of the Feynman- 
Vernon influence functional technique^ Then the parti- 
tion function is reduced to a single functional integral 
over x, i.e., Z = j Vx(t) exp (S 2 g[x]/h), where the 

effective action is given by S^§[x] 
with 36 i 37 
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dr I dr' [x(t) - x(r')f K(t - r') 



(11) 

Thus the dynamics of the phase difference with an 
anisotropic mass moving in a two-dimensional potential 
U(x,y) can be mapped into a simple one-dimensional 
model. Note that due to the coupling between the junc- 
tion and the LC circuit, the effective action S 2 $ [x] con- 
tains a kind of dissipation action [x] in a sense that 
a retardation (or nonlocal) effect exists. 

The nonlocal kernel K(t) in Eq. (fTTj) is defined by 
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where wlc = 1/V-^sCs is the LC resonant frequency, 
to n = 2im/hf3 is the Matsubara frequency, and f„ = 
(Cs/Cj)uf jC /(uj 2 + uj\ c ) is the Fourier coefficient of 
K(t). The nonlocal kernel K(r) is related to the spec- 
tral density J (to) in the language of the Caldeira-Leggett 
theory 3 ^ via 
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where -D w (t) is the Matsubara Green's function of 
a free boson, £> w (r) = (l/ft/3) J2™=-oo GX P («w„r) 
{2to/(to 2 + to 2 )} . In the Caldeira-Leggett theory, the ex- 
ternal bath is modeled by a sum of an infinite num- 
ber of harmonic oscillators which is capable of destroying 
MQT. 3,3 - On the other hand, in our case, the bath can 



be described by a single harmonic oscillator and then 
J (to) is given as a delta-function, i.e., 



J(w) = -mto\ c 6(u) - to LC ). 



(14) 



For later convenience, we rewrite the effective action 
S*g[x] in terms of the memory kernel k(r),— i.e., 
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where the Fourier coefficient of the memory kernel is 
related to the Fourier transform of the memory-friction 
kernel j(t)£ i.e., j(to) = J™ dtj{t)e' A 



via 
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In below, by using the derived effective action S 2 ^ , we 
calculate the crossover temperature and the escape rate, 
and then discuss the influence of the anisotropic mass 
and the two-dimensional nature of the potential profile on 
the macroscopic phase dynamics in the YBCO junction 
coupled to the LC circuit. 



IV. CROSSOVER TEMPERATURE 

Based on the normal mode analysis at the barrier top of 
the potential U (x, y), the crossover temperature between 
the MQT and TA escape process is given by 3 ^ 
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where ujr is the renormalized trial frequency due to the 
coupling to the LC circuit and is a positive root of the 
equation, + wr7(wr) = to 2 . By using Eq. (jTTJ) , we get 




for the adiabatic ( Lj ^S> L$ and Cj Cs ) and the 
nonadiabatic (Lj -C L$ and Cj <C Cs ) limit, where 
9 = 1- (l + C s /Cj)to 2 c /to 2 p , and Lj = Ljo/^T^ 2 
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is the Josephson inductance. In the adiabatic (nonadi- 
abatic) limit, the LC resonant frequency c^lc is much 
larger (smaller) than the Josephson plasma frequency u> p . 
As a result, the LC resonant mode can (cannot) adiabat- 
ically follow the dynamics of the Josephson junctions. 
Note that, in the YBCO biepitaxial junction used in the 
MQT experiment the nonadiabatic case is approxi- 
mately realized (see Sec. VII). 

By substituting Eq. (|20|) into Eq. fTS]) , we get a simple 
expression for T co as 
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where = Huj p /2 irks is the crossover temperature 

without coupling to the LC circuit. Therefore, due to 
the influence of the LC circuit, T^P becomes smaller than 
Tjjp. In the nonadiabatic (adiabatic) case, reduces 
with decreasing the kinetic inductance Ls (increasing the 
stray capacitance Cs). 



THERMAL ACTIVATION PROCESS 



The TA escape rate well above the crossover tempera- 
ture T^P is given by^ 



t^2D 
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where CqP is the quantum-mechanical enhancement fac- 
tor resulting form stable fluctuation modes and j(u) n ) = 
(Cs/Cj)\uin\u\ c /{LOn + oj\ c ) is the Fourier transform of 
the memory-friction kernel. As was shown in Sec. II, 
the potential barrier height is not changed even in the 
presence of the LC circuit, i.e., V 2D = V 1D . Therefore 
the coupling to the LC circuit only modifies the prefactor 
of Tta- Note that the exponent in Tta for the Joseph- 
son junction coupled to the LC circuit was discussed by 
FistuL 40 Here we calculate the TA escape rate r 2 ^ in- 
cluding the prefactor in order to discuss the influence of 
the LC circuit explicitly and compare with experimental 
results. The quantum correction c qm in the prefactor can 
be calculated analytically as 
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where p± = (1 + C s /Cj) (w L c/wi) 2 ±(w p /wi) 2 , and q = 
WpWLc/^i. The effective trial frequency in Eq. I|22p 
is given by Eq. (fl9|) . 

In the adiabatic limit (Lj S> Ls and Cj 3> Cs), the 
quantum enhancement factor c^P (|24p can be simplified 
to 
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On the other hand, in the non-adiabatic limit (Lj < Lj 

and Cj < Cs), we get 
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Note that the nonadiabatic limit is almost realized in 
the actual YBCO junction^ (see Sec. VII). The quan- 
tum enhancement factor c 2 P almost coincides with the 
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result without retardation effects^ i.e., c^ m 
smh(h/3u) p /2) I sm(h/3ujp/2) for both the adiabatic and 
non-adiabatic limit (see also Sees. VI. C and VII. B). 
From Eqs. fl2UJ), ([27]), and J2SJ, therefore, we can 

conclude that the influence of the coupling to the LC 
circuit on the thermal activation process is quite weak, 
so the system behaves as a one-dimensional system well 
above the crossover temperature . This result is qual- 
itatively consistent with the experimental observation in 
the YBCO biepitaxial junction^ In Sec. VII we will nu- 
merically compare theoretical results with experimental 
data in the TA regime. 



VI. MACROSCOPIC QUANTUM TUNNELING 
PROCESS 



To obtain the MQT escape rate T^rp for the Joseph- 



son junction coupled to the LC circuit below T C( P, the 
usual procedure is to apply the so-called Im-F mcthodi^ 
It is based on the calculation of the free energy F = 

— (l//3)lnZ and thus of the partition function Z of an 
unstable system. According to the metastable decay 
theory^ the MQT escape rate has the form T^-p = 

— (2/ft)ImL. In below we derive a weak retardation 
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in the effective action S^g can 
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be treated as a small perturbation. Then we calcu- 
late r^Q T for zero temperature and finite temperature 
(0 < T < T C 2 D ) in a weak retardation limit. 



A. Weak retardation condition 

Here we derive a condition for the weak retardation 
or nonlocal effect. The bounce trajectory xb(t) satisfies 
the Euler-Lagrangc equation obtained from the variation 
principle 6S%g[x] = 0. Taking the variation of the effec- 
tive action Eqs. (flO|) and (fTTj) . we then get 



M 



d 2 x B (t) _ dUxp 
dr 2 8xb(t) 



Mi 

\ dT , k(T-T , )x B (T , ) = 0.(29) 
Jo 



Introducing the Fourier expansion xb{t) = Y^=-oo Xn 
exp(ia; n r) into Eq. (j2"5)) and using Eq. p^)l . we obtain 
the following equation for the Fourier coefficient i„ 
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For uj n = Up, we get 
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The criterion for weak retardation can be obtained from 
Eq. (|3~Tj) by comparing the retardation term (M£„ (u>„ = 

Wp)iL B ' 1 ) with the first term (2Mw 2 i:„ ). By defining 
the parameter— 
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Cj < C s 



the weak retardation limit corresponds to the case of 
k <C 1. In this case 5^1* in the effective action can be 
treated as a small perturbation. Figure 3 shows the Ls 
and Cs dependence of the retardation parameter n. In 
the YBCO junction used in the MQT experiment r*^- 
Cs/Cj ~ 7.3 and Lg/Lj w 3.2 as will be estimated in 
Sec. VII, so n w 0.15. Therefore, the weak retardation 
condition is almost satisfied in this case. 
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FIG. 3: (Color online) Ls and Cs dependence of the pa- 
rameter k which characterizes the retardation effect due to 
the junction-LC coupling. In the YBCO biepitaxial junction 
which was used in actual MQT experiment, n « 0.15 (red 
circle) . 



techniques, the MQT escape rate r^Q T in the weak re- 
tardation limit k <C 1 is perturbatively determined by2£ 



rM QT ( T = °) = -Vl20^ (Bid + B lct ) e - B ^- B ^ ,(33) 
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where Bid = S^g [xb]/7i = 36Vo 1D /5fta;p and B lct = 
Slg [x B ]/fi arc the bounce exponents, that are the value 
of the actions evaluated along the bounce trajectory 
xb(t) = x isech 2 '(u> [ p t / '2). The bounce action S^[xb] 
can be written as 
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B. Zero Temperature 



Therefore the total bounce exponent is given by 



The MQT escape rate at zero temperature is given by 
lim^^ 00 (2//3)ImlnZ4 By use of the bounce 



p2D 
1 MQT 



B\T> + -Bdiss 



15ft 



(M + 5M)w p x\, 



(36) 



7 



where 
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L s Jo 



z 2 + ^ } sinh z (ttz) 



5 Lj > ( 37 ) 

is the retardation correction to the mass M. Thus, due 
to the influence of the LC circuit, the bounce exponent 
is increased with respect to the one-dimensional case. 

By substituting Eq. (j3l)j) into Eq. (j3"3")l . we finally get 
the zcro-tcmpcraturc MQT escape rate for c < 1 as 
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Therefore the coupling to the LC circuit effectively in- 
creases the barrier height, i.e., V 1B -> V 1D (1 + SM/M). 
This behavior is consistent with the result for low-T c 
junctions coupled to an LC circuit which was derived in 
a different contexti 37 ' 42 i 43 In contrast to the TA regime, 
the coupling to the LC circuit reduces the MQT escape 
rate L^q T considerably. Therefore the anisotropic mass 
and the two-dimensional nature of the potential profile 
have large influence on the MQT escape process in the 
low-temperature regime. 



C. Finite Temperature 

In this section we calculate the MQT escape rate at 
finite temperature (0 < T < T£?) for the weak retar- 
dation limit (k « 1) and show that the T dependence 
of the MQT escape is drastically influenced by the cou- 
pling to the LC circuit. At finite-temperature the MQT 
escape rate is enhanced by thermal fluctuations. If the 
temperature is much lower than Tlc = hu>Lc/2irkB, the 
finite temperature enhancement to the r^g T (T = 0) is 
negligible in contrast to the system with an Ohmic dissi- 
pative environment (sec Appendix A). As will be shown 
later, however, at sufficiently high temperature but still 
smaller than T C 2 Q D , the coupling to the LC circuit gives a 
large enhancement to the r 2 ^Q T (T = 0). 

The finite-temperature bounce trajectory x^(t) is 
given by a periodic solution in the inverted potential 
-U m {x) with energy -E (0 < E < V 1D ). The cor- 
responding solution is given by x%(r)/x\ = q2 + (qi — 
qz) cn 2 (Aw p T|fc), where < q2 < qi are the solutions 
of the equation q 2 (l - q) = AE/27V m M. In this equa- 
tion cn denotes the Jacobi elliptic function with the pa- 
rameter k = y(<7i — 52)/ {qi — 93) and the coefficient 
^ = V c h ~ 93/2. The period h(3 of the finite-temperature 
bounce x^ (t) is given by the complete elliptic integral of 



the first kind K(k) = J* dx [(1 - x 2 )(l - kx 2 )} 1/2 as 
hf3 = 4,K{k)/u} p ^Jqi — q 2 ^ The energy E is related to 
the temperature /3 _1 via this equation. In the limit of 
zero-temperature (E — > 0), we recover the zero tempera- 
ture bounce solution xb(t) = xisech 2 (w J ,T/2). 

By substituting i| (r) into the effective action, a finite- 
temperature bounce action iS|P[:eb] can be evaluated as 



S^P[x B }+S r T ct [x B 
M 

lo 

\uJ p M{q 1 -q 2 ) 2 



(39) 



ftp 

/ dTx T B (rf + Ehf3 
Jo 

F(Xuj p hf3, k) + EH/3, (40) 



15fc 2 



1 r h0 

SF\x B ] = g J o dTdr'K(T-T')[x T B {T)~x T B {r')\ 



= —m 



LC 



n— — 00 



LO 2 + LU 



\x T B {n)\ , (41) 



LC 



where x^ (n 



(l/hf3)f Q drx^ (t) exp (iw„t) is the 
Fourier transform of the finite temperature bounce x 1 ^ (r) 
and the function F(c, k) is defined by 



F(c, k) 



dn(c|fc) 



y/1 - fcsn(c|fc) 2 



[d(k) + e(k)E(am(c\k),k) 



fcsn(c|fc) cn(c|fc)/(c, k)\Jg(c, k) , 



(42) 



with sn(c|fc) and dn(c|/c) the Jacobi elliptic func- 
tions, am(c|fc) is the Jacobi amplitude, E{4>,k) = 

Jq d9y/ 1 — k 2 sin 2 9 is the elliptic integral of the second 
kind, d(k) = -4 {2 + k(k - 3)}, e(k) = 8 {1 + k(k - 1)}, 
f(c,k) = k-2-3kcn 2 {c\k) + 3ksn 2 (c\k), and g(c,k) = 

2kcn 2 {c\k) - 2 {k - 2 + k sn 2 (c|fc)} 2 . 

The prefactor of the MQT escape rate r 2 ]^Q T (T) is de- 
termined from the fluctuation modes around the bounce 
solution x%(t) and can be calculated from the Gel'fand- 
Yaglom formul: 



,45,46 



A(T) 



1 

2 

sinh 



9M W 3(1 - 7 2) {qi _ q 3 )V^ qi - q 2 )(l - k 2 ) 



2irh 

u>p/3h 



y/a{k)E{k) 



b(k)K(k) 

(43) 



where a{k) = 2(fc 4 — k 2 + 1), b(k) = (1 - k 2 ){k 2 - 2), and 
E(k) is the elliptic integral of the first kind41 For E — * 0, 
A(T) is reduced to well-known zero-temperature result, 
i.e., A(T = 0) = (uj p /2n) v /864:TTV lu /hLjp. Thus we ana- 
lytically obtain the finite-temperature MQT escaper rate 
for the weak retardation limit: 



^2D 
MQT 



(T) = A(T)exp 



S™{x B ] 



(44) 



This result is smoothly matched to the zero-temperature 
expression ([3"5t when T — > 0. 
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The Arrhenius plot of the escape rate T|[q T (T) and 
r^(T) are shown in Fig. 4. Note that we have used 
an interpolate crossover-formula 3 " near the crossover tem- 
perature T C 2 Q D because Tff QT {T) and T^(T) diver ges as 
T — > T C 2 D D due to the breakdown of the Gaussian approx- 
imation (see Appendix B). In the numerical calculation, 
we have used Ic = 1.4 /iA and Cj = 0.22 pF. These 
parameters agree with those directly measured or esti- 
mated from experiments^ allowing a comparison be- 
tween numerical and experimental data (see Sec. VII). 
The TA escape rate r|^(T) above T 2 D D is almost same 
as in the case without the LC circuit as clearly seen from 
Fig. 4. On the other hand, r 2 ^Q T (T) is considerably 
reduced due to the coupling to the LC circuit. As was 
explained in Sec. II, the MQT escape rate is significantly 
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FIG. 4: (Color online) Arrhenius plot of the escape rate L(T') 
for the YBCO junction coupled to the LC circuit, (a) Ls 
dependence of T(T) for Cs = L6 pF. (b) Cs dependence of 
r(T) for L s = 1.7 nH. We choose 7 = 0.9, I c = 1.4 fxA (Lj = 
0.54 nH) and Cj = 0.22 pF for both graphs. In the above 
parameter range, the weak retardation condition (k <C 1) 
is satisfied. Solid and dotted arrows indicate the crossover 
temperature T C 2 D and Tlc, respectively. 



reduced with decreasing L5. This is due to the drastic 
change in the two-dimensional potential profile. Simi- 
larly the MQT escape rate is decreasing with increasing 
Cs- This behavior can be explained as follows. In the 
case of small Cs, the kinetic energy to the y (or 0s) 
direction becomes small. Therefore, due to the strong 
anisotropy of the mass, MQT along the external escape 
direction which connects the minimum (4> m , 0™) and the 
saddle point ((/>*, </> s ) (see Fig. 2) is inhibited. 

Interestingly, the temperature dependence of the MQT 
escape rate r 2 ^Q T (T) for the two-dimensional system is 
quite different from that of the one-dimensional system. 
In the case without a LC circuit, FJJq T (T) is almost con- 
stant below the crossover temperature . On the other 
hand, in the case of small Ls and large Cs, rMQT^) 
significantly depends on the temperature in the range of 
Tlc < T < T" C 2 D . Therefore, from the temperature de- 
pendence of the MQT escape rate, we can distinguish 
qualitatively the anisotropy of the mass and the dimen- 
sionality of the quantum phase dynamics. 



VII. COMPARISON WITH EXPERIMENTS 

In this section, we first summarize the experimental de- 
vice parameters relevant to the TA and MQT phenomena 
of a YBCO grain-boundary biepitaxial Josephson junc- 
tion (Sec. VII. A). Then, in order to check the validity of 
the extended circuit modet!^ 3 .^ for the TA and MQT 
escape process, we try to compare our result with the ex- 
perimental data£ of the switching current distribution at 
the high-temperature TA and the low temperature MQT 
regimes (Sec. VII. B). 



A. Experimental parameters 

Experimental parameters used in numerical calculation 
are given in Table I. The values given in the table are in- 
deed typical ones in actual experiments. By the fitting 
to the Kramers formula with experimental data of the 
TA escape rate well above the crossover temperature, Ic 
was experimentally estimated as Ic = 1.4 /iA which cor- 
responds to Ljo = 0.24 nH£ The values of Ls = 1.7 nH 
and Cs = 1.6 pF have been directly determined from 



Parameters 


Estimated value 


Ic- Josephson critical current 


1.4 rtA 


Ljo: Josephson inductance 


0.24 nH 


Lj(-y = 0.9): Josephson inductance 


0.54 nH 


Cj: Capacitance of the junction 


0.22 pF 


Ls: Inductance of the LC circuit 


1.7 nH 


Cs'- Capacitance of the LC circuit 


1.6 pF 



TABLE I: Experimental junction parameters used for the es- 
timation of the TA and MQT escape rate. 
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the ELQ experiments^ Remaining unknown junction 
parameter is Cj. As will be shown in the next section, 
this value is estimated as Cj — 0.22 pF. This allows for 
a qualitative comparison between numerical and experi- 
mental data of the zero temperature MQT escape rate. In 
this case Cs/Cj w 7.3 and L 5 /Ljo « 7.2 (L s /Lj « 3.2 
for 7 = 0.9). Therefore, in the actual junction, the nona- 
diabatic limit is approximately realized. 



B. Numerical results 

In this section, we numerically calculate the switching 
current distribution P(j) which is related to the escape 
rate T as^- 



P( 7 ) = ~r( 7 )exp 



r( 7 W 



(45) 



where v = \dn/dt\ is the sweep rate of the external bias 
current. In the actual experiment,— the temperature de- 
pendence of the full width at half maximum (HMFW) a 
of P(j) is measured as shown in Fig. 5. 

First we investigate the TA regime. In Fig. 5 we 
show the temperature dependence of a in the TA es- 



cape regime (red solid and black dotted lines). In calcu- 
lation we have substituted Eqs. (|22|) and (|24|) into Eq. 
(|45| . Both the one- and two-dimensional model give good 
agreement with the experimental data (black circles) well 
above the crossover temperature (T C 2 D D ~ 0.09K). There- 
fore, in the TA regime, the system can be treated as a 
one-dimensional model without the LC circuit. 

In the MQT regime, the measured saturated value of 
cr at T = 0.03 K is found to be 11.9 nA£ From the 
numerical estimation of cr, we found that Cj — 0.22 pF 
gives good agreement with the experimental value of cr as 
shown in Fig. 5. The obtained value of Cj is consistent 
with the estimated value Cj « 0.16 pFsa O.lCs based 
on the geometry of the junction . 33 ' 50 Therefore, we can 
conclude that the extended circuit model quantitatively 
explain the MQT experiment 5 - in the YBCO biepitaxial 
junction near zero temperature. 

In order to test the validity of the extended circuit 
model more systematically, experimental measurements 
of HMFW a by changing L$ and Cs are needed.— In 
Fig. 6 we show L$ and C's dependence of cr. If we use 
the substrates with low dielectric constant e, e.g., MgO 
(e ~ 9.6) and LaAlC>3 (e ~ 23), the reduction in a with 
respect to the one-dimensional model (blue dotted line 
in Fig. 6) becomes small. 
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FIG. 5: (Color online) The temperature dependence of the 
full width at half maximum a of the switching current dis- 
tribution P(7). Both the one- (dotted black line) and two- 
dimensional models (red solid line) give almost same result 
above the crossover temperature. The calculated a from the 
zero-temperature MQT escape rate for the one- dimensional 
model and the finite-temperature MQT escape rate for two- 
dimensional model [Eq. (44)] are shown by dashed-dotted 
(blue) line and (red) squares, respectively. Experimental data 
of a (black circles) for a YBCO biepitaxial Josephson junc- 
tion^ is also plotted. We choose la = 1-4 fiA (Ljo = 0.24 
nH), Cj = 0.22 pF, L s = 1.7 nH, C s = 1.6 pF and vl c = 1.0 
mA/s. Arrows indicate the several characteristic tempera- 
tures, i.e., T LC , T™, and T C 2 D for 7 = 0.9. 



VIII. CONCLUSIONS 

In the present work, the TA and the MQT escape 
process of the YBCO Josephson junction coupled to the 
LC circuit has been analyzed by taking into account the 
anisotropy of the mass and the two-dimensional nature 
of the phase dynamics. Based on the Feynman- Vernon 
approach, the effective one-dimensional action is derived 
by integrating out the degree of freedom of the LC cir- 
cuit. We found that the coupling to the LC circuit gives 
negligible reduction for TA escape rate. On the other 
hand, we also found that the MQT escape rate is consid- 
erably reduced due to the coupling between the junction 
and the LC circuit. More importantly, the temperature 
dependence of the MQT escape rate for the YBCO junc- 
tion coupled with the LC circuit is quite different from 
that without the LC circuit. These theoretical results 
are in an excellent agreement with experimental data of 
the YBCO biepitaxial Josephson junction* 5 - Therefore we 
can conclude that the anisotropy of the mass and the 
two-dimensional nature of the potential profile due to the 
coupling to the LC circuit are quite important and essen- 
tial to understand macroscopic quantum phenomena and 
qubit operation in such systems. 

There, however, remain a question, which cannot be 
treated by the present approach. As seen from Fig. 5, the 
extended circuit model cannot explain a hump structure 
of a near 0.1 Kj 5 -^ We note that such a characteristic 
behavior has been observed in a DC-SQUID (supercon- 
ducting quantum interference device) system composed 
only by low-T c superconductors £1 While the deviation 
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FIG. 6: (Color online) (a) L s and (b) C 3 dependence of the full 
width at half maximum a of the switching current distribution 
P(7) at T — K for the case with (red solid) and without 
(blue dotted) the LC circuit. We choose the parameters as 
I c = 1.4 fiA (Lj = 0.24 nH), Cj = 0.22 pF and vie = 1-0 
mA/s. 



per can be regarded as an artificial atom (the Joseph- 
son junction) coupled to the quantized electromagnetic 
field (the LC circuit). Therefore the appearance of sev- 
eral interesting phenomena relating to quantum optics^ 
e.g., the vacuum Rabi oscillation^ generation of a non- 
classical state of the LC system^ and laser oscillation^ 
is expected also in the YBCO biepitaxial junctions. Ad- 
ditionally the LC circuit will act as a quantum informa- 
tion busi 58 i 59 Therefore the entanglement or the coupling 
between separated high-T c qubits and eventually a high- 
Tc version of the circuit-QED syste m 60 i 61 ' 62 will be re- 
alized in such biepitaxial junctions. These studies will 
open up the possibility of future applications for high-T c 
superconductor materials. 
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APPENDIX A: FINITE TEMPERATURE 
CORRECTION OF THE MQT ESCAPE RATE 



between our theory and the experimental data is left for a 
future problem, it might be attributed to the fact that our 
model neglects thermal/dynamical population from the 
quasi-ground state to the excited states in the metastable 
well. We note that in the nonadiabatic cases (ui p ^> i^lc) 
corresponding to the MQT experiment for the YBCO 
biepitaxial junction^^ more quantum levels due to ex- 
citation of the LC circuit are relevant to the decay pro- 
cess than in a simple one-dimensional system. In order 
to investigate the thermal/dynamical population effect, 
we have to solve a master equation to obtain population 
probabilities of each level by the Larkin and Ovchinnikov 
theory^ This consideration may explain the anomaly of 
the escape rate T(T) near the crossover temperature. 

Finally we would like to comment advantages of the 
junction-LC coupling to the qubit and quantum optics 
applications. The system that we considered in this pa- 



If the temperature T is much smaller than the crossover 
temperature T™, the finite-temperature MQT escape 
rate (in the weak retardation limit) is given by the prod- 
uct of r^Q T (T = 0) and a finite-temperature correction 



.,63 



MQT 



(T) 



^2D 
MQT 



(T = 0)exp[A(T)}, (Al) 



where A(T) is the finite-temperature correction to the 
bounce action and is given in term of the spectral density 
J(w) as 



A(T) 



2irh 



dusJ((jj) 



coth 



h(3uj 



(A2) 



In this equation xbtb = /^(iraBtr), where x_b(t) is 
the zero temperature bounce. For the spectral density 
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Eq. (14), the enhancement function takes the form 
A{T) 



I i 



!^|0 (1 _ 7 2 )e -^ LC; 



(A3) 



for T <C fvjJLc /2nkB = Tlc- Thus the thermal enhance- 
ment shows exponentially weak temperature dependence 
as long as T <C TLc- On the other hand, for a damped 
system with J (to) oc lo s , we get algebraic large enhance- 
ment A(T) oc T 1+S M&$. Thus the coupling to the LC 
circuit gives the weak finite-temperature correction to 
the zero-temperature MQT escape rate r^Q T (T = 0) 
for T < Tlc- Note that, in the actual YBCO junc- 
tion, T LC < T C 2 D . In Sec. VI. C, we show the finite- 
temperature MQT escape rate rfJq T (T) for < T < 

rp2D 



APPENDIX B: THE ESCAPE RATE IN THE 
CROSSOVER REGIME 
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FIG. 7: (Color online) Temperature dependence of the escape 
rate near the crossover temperature T C 2 D . The TA escape rate 
TTA(r) is represented by a dotted (red) line and the MQT 
escape rate r|Jq T (T) as a dashed (blue) line. The crossover 
formula r 2 ^ (T) (B3) is represented by a solid (black) line and 
smoothly matches onto these functions below and above the 
crossover temperature T C 2 D (arrow). Parameters are Ic = 1.4 



In this Appendix, we show an expression for the escape ^a, Cj = 0.22 pF, L s =1.7 nH, and C s = 1.6 pF 



rate near T 2 ^ . The quantum enhancement factor c 2 ]^ in 
the TA escape rate T^(T) increases with decreasing the 
temperature, and diverges as T -> T C 2 G D . This unphysi- 
cal divergence is due to the appearance of a bounce tra- 
jectory for temperature below T 2 ^ and hence the break- 
down of the Gaussian approximation^ A crossover region 
is characterized by the condition |T — T C 2 D | < T 2 D /<5 in 
which the dimensionless parameter 8 is given by^ 



1 



di{uj R ) \~\ [MP 
du R )\ V 2B 4 



(Bl) 



Here the coefficient B4 measures the strength of the an- 
harmonicity of the potential Uid and is given by 



Ba = uj. 



pO 



1 



1 



u p0 



r 



2w| 



w 2 +^27(^2) 



.(B2) 



In the crossover region, functional integral cannot be 
done by steepest descents, but requires a more careful 
treatment, leading to a escape rate of the form^ 



is the dimensionless prefactor, and |r'| = H/3 2 /2k(P = 
/? co ) 2 is the energy derivative of the bounce period. 
This expression applies for the temperature T slightly 
above and below T 2 ^ and smoothly interpolates between 
the high-temperature escape rate r 2 J^(T) and the low- 
temperature escape rate T 2 |q T (T) as demonstrated in 
Fig. 7. 
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with erfc(x) = exp(— y 2 /2)/\/27r. In this equa- 

tion, /3 co = 1/fcsT- 
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